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Abstract 

We allow representing and reasoning in the 
presence of nested multiple aggregates over 
multiple variables and nested multiple aggre- 
gates over functions involving multiple vari- 
ables in answer sets, precisely, in answer set 
optimization programming and in answer set 
programming. We show the applicability of 
the answer set optimization programming with 
nested multiple aggregates and the answer set 
programming with nested multiple aggregates 
to the Probabilistic Traveling Salesman Prob- 
lem, a fundamental a priori optimization prob- 
lem in Operation Research. 



1 Introduction 

A fruitful approach to aggregates in answer set program- 
ming has been presented in |Faber et al., 2010] that al- 
lows to represent and reason in the presence of aggregate 
function which are defined through aggregates atoms 
that are allowed to be recursive. The aggregate functions 
defined in |Faber et al., 2010] , and hence the aggregate 
atoms defined over them, are allowed to be monotone, 
anti-monotone and non-monotone which is clearly ap- 
propriate for non-monotonic reasoning. The answer set 



semantics with aggregates defined in Faber et al, 2010 



have shown to be appropriate for many interesting prob- 
lems in answer set programming that require aggregates 
over some numerical criteria imposed by the problems 
and arise in many domains. 

On the other hand aggregates have been consid- 
ered in answer set optimization programming described 
in |Saad and Brewka, 201 1| that extended the answer 
set optimization programs |Brewka et al., 2003] to allow 
representing and declaratively solving multi-objective 
optimization problems in an answer set programming 
framework. The answer set optimization with aggregates 
framework of |Saad and Brewka, 2011 have shown to be 
able to find Nash equilibrium in strategic games with any 
number of players and with any number of strategies, 
which is a multi-objective optimization problem with 
multiple conflicting goals. 



However the aggregates defined in the answer set pro- 
gramming framework in |Faber et al., 2010] allow a sin- 
gle aggregation over a single variable which limits its ap- 
plicability to many interesting problems, especially the 
problems that require aggregations over a function with 
multiple variables. And hence, more generally limits its 
applicability to the problems that require nested mul- 
tiple aggregations over multiple variables and the prob- 
lems that require nested multiple aggregations over func- 
tions with multiple variables. 

On the other hand, although the aggregates 
defined in answer set optimization programming 
|Saad and Brewka, 2011] , is similar to the aggregates de- 
fined in the answer set programming [Fab er et al, 2010 , 
allows a single aggregation over a single vari- 
able, the answer set optimization programming in 
|Saad and Brewka, 2011] is still capable of solving some 
interesting multi-objective optimization problems like 
the one arises from finding Nash equilibrium in strate- 
gic games. However, the answer set optimization pro- 
gramming of | Saad and Brewka, 2011 is still incapable 
of representing and reasoning about many interesting op- 
timization problems that require optimization over ob- 
jective functions that involve nested multiple aggrega- 
tion over functions that contain multiple variables. This 
implies minimization or maximization over nested mul- 
tiple aggregations over functions with multiple variables 
or nested multiple aggregations over multiple variables. 

In this paper we generalize both answer set optimiza- 
tion programming |Saad and Brewk a" 2011] and answer 
set programming [Faber et al., 2010 with nested mul- 
tiple aggregates to allow representing and reasoning in 
the presence of nested multiple aggregation over multi- 
ple variables and nested multiple aggregation over func- 
tions with multiple variables. We show the applicabil- 
ity of the nested multiple aggregates in answer sets to 
the Probabilistic Traveling Salesman Problem (PTSP), a 
fundamental stochastic optimization problem in Oper- 
ation Research |Jaillet, 1988] , whose objective function 
involving nested multiple aggregations over a function 
with multiple variables. 

Probabilistic Traveling Salesman Problem is a priori 
optimization problem for which we want to find a priori a 
tour with minimum (expected) length through a set of n 



points, where only k points out of n points (0 < k < n) 
at any given instance of the problem must be visited, 
where the number k is chosen at random with a known 
probability distribution, that is based on the probability 
distribution of visiting each point in the set of n points. 
In addition, the k points chosen to being visited, in a 
given instance of the problem, has to be visited in the 
same order as they appear in the a priori tour. 

Formally, the Probabilistic Traveling Salesman Prob- 
lem is defined as follows. A Probabilistic Traveling 
Salesman Problem (PTSP), G, is a tuple of the form 
G = {V, E, D, P), where (V, E) is a graph with a set of 
n nodes, V, and a set of arcs, E, connecting the nodes 
in V , D is an n x n (complete) matrix (the distance ma- 
trix) that specifies the distance (cost), dy, incurred by 
traveling from a node (location) i to a node (location) 
j, P = {pi,p2, ■ ■ ■ ,p n } is a probability distribution over 
the nodes in V, where pi is the probability of a node i 
must be visited. In other words, the nodes in V repre- 
sent the set of locations that need to be visited, arcs in E 
represent the roads connecting the locations, the values 
dij in D represent the distance traveling from location 
i to location j, and the probabilities pi in P represent 
the probability that a location i must be visited. The 
objective is to find a priori tour t, a Hamiltonian cir- 
cuit of G, with the minimum expected length (cost). By 
indexing the nodes of a tour t by their order of appear- 
ance, we present the tour, t, as the sequence of nodes 
t = (0, 1, 2, 3,4, . . . , n,0), where is the depot or the 
starting point from which the tour always starts and 
ends. It has the index and with probability pa = 1, 
since it must be always visited. Therefore, the expected 
length, L t of a tour t, denoted by E[L t ], is given by: 
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where i,j, k are indices of nodes within a tout t (not the 
nodes themselves) and n+l = 0. Given that T is the 
set of all tours in G, the objective is to find a tour t in 
T with the minimum expected length E[L t ], i.e., solve 
the a priori optimization problem 

min E\L t ]. 
teT 

In this paper we develop syntax and semantics for both 
answer set optimization programs and answer set pro- 
grams to allow representing and reasoning in the pres- 
ence of nested multiple aggregates over multiple vari- 
ables and nested multiple aggregates over functions in- 
volving multiple variables. We show the applicabil- 
ity of both the answer set optimization programs with 
nested multiple aggregates and the answer set pro- 
grams with nested multiple aggregates to the funda- 
mental stochastic optimization Probabilistic Traveling 
Salesman Problem, a priori optimization problem in Op- 
eration Research Ijaillet, 1988 . In addition, we show 



optimization programs with nested multiple aggregates 
and in answer set programs with nested multiple ag- 
gregates, a corresponding classical traveling salesman 
problems can be intuitively represented and solved in 
both frameworks. We prove that the presented an- 
swer set optimization programs with nested multiple 
aggregates modify and generalize the answer set opti- 
mization programs with a single aggregate over a sin- 
gle variable proposed in [Saad and Brewka, 201 1| as well 
as a generalization of the answer set optimization pro- 
grams described in |Brewka et al, 20 03 . In addition, 
we prove that the presented answer set programs with 
nested multiple aggregates generalize the answer set pro- 
grams with a single aggregate over a single variable de- 
scribed in |Faber et al., 2010] as well as a generaliza- 
tion of the original answer set programs presented in 
|Gelfond and Lifschitz, 199T] . 

2 Nested Aggregates Disjunctive Logic 
Programs 

In this section we present the syntax and the answer 
set semantics of disjunctive logic programs (a form of 
answer set programming) with nested multiple aggre- 
gates, denoted by NDLP. The syntax and the answer set 
semantics of NDLP generalize and modifies the syntax 
and semantics of disjunctive logic programs with a sin- 
gle aggregate over a single variable, DLP- 4 , presented in 
|Faber et al, 2010 . 



that with minor modification to the probabilistic trav- 
eling salesman problems representation in answer set 



2.1 NDLP Programs Syntax 

Let £ be a first-order language with finitely many predi- 
cate symbols, function symbols, constants, and infinitely 
many variables. A term is a constant, a variable or a 
function. An atom is a predicate in Be, where Be is the 
Herbrand base of C. The Herbrand universe of C is de- 
noted by Uc- Non-monotonic negation or the negation 
as failure is denoted by not. 

An expression of the form {F | G} is called a sym- 
bolic set, where F is a single variable or a function F — 
g(Xi, . . . , X n ), where X\, . . . , X n are variables or func- 
tions and G is a conjunction of atoms, negation of atoms, 
aggregate atoms, and negation of aggregate atoms (de- 
fined below). However, a ground set is a set of pairs of 
the form (Fg | Cg) such that Fq is a constant and Cg is 
a conjunction of ground atoms, negated ground atoms, 
ground aggregate atoms, and negated ground aggregate 
atoms (defined below). A ground set or a symbolic set is 
called a set term. We say f(S) is an aggregate function 
if / is an aggregate function symbol and S is a set term, 
where / G {min, max, count, sum, times}. If f(S) is an 
aggregate function and T is a constant, a variable or a 
function term called guard, then we say f(S) -< T is an 
aggregate atom, where -<G {=, +S <, >, <, >}. 

Observe that, from the definition of aggregate func- 
tions and aggregate atoms, it can be seen that an aggre- 
gate function f(S) in an aggregate atom f(S) ~< T can 
be defined in terms of other aggregate atoms f'(S') -< T" 
and hence in terms of other aggregate functions f'(S'). 



This means multiple level of nested aggregate functions 
calculation is achieved by allowing aggregate functions 
of the aggregate atoms to be defined in terms of other 
aggregate atoms. This is accomplished by allowing the 
conjunction in the definition of the symbolic set S of the 
aggregate function f(S) in the aggregate atom f(S) -< T 
to contain other aggregate atoms as constituents. This 
can be illustrated by the following example. 

Example 1 Assume that we want to represent the dou- 
ble summation 



J2J2(*+j) 



(2) 
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as an aggregate function, f(S), and assign the result to a 
variable called X to form the aggregate atom f(S) = X . 
Assume also that the possible values of i,j are repre- 
sented by the predicate d{I, J) . Therefore, the double 
summation (0) can be represented as the aggregate atom 
f(S) = X that is defined as: 

sum { A | 1 < I < n, 
sum { 1+ J | d(I, J), 1 < J < m, } = A } = X (3) 

where S={A\l<I<n, sum { I + J \ d(I, J), 1 < 
J <m, } = A}. 

Notice that there are two instances of the aggregate func- 
tion sum in the aggregate atom representation of the 
double summation ^ . The outer instance of the aggre- 
gate function sum and the inner instance of the aggre- 
gate function sum. The inner instance of the aggregate 
function sum is defined inside the outer instance of the 
aggregate function sum. 

Moreover, notice that the variable J occurs only in the 
inner instance of the aggregate function sum. Whereas 
the variable / occurs in both the outer and the inner in- 
stances of the aggregate function sum. This means that 
the variable J is seen only by the inner aggregate func- 
tion sum, but not seen by the outer aggregate function 
sum, and I is seen by both the outer and the inner ag- 
gregate function sum. This implies that the occurrence 
of I in the outer sum hides the occurrence of / in the 
inner sum, which means that any substitution of / to a 
constant in both the inner and the outer sum is made 
by the outer sum. This also implies that the scope of 
the variable J is the inner sum while the scope of the 
variable / is both the outer and the inner sum, i.e., the 
entire symbolic set of the outer sum. This motivates the 
following definition of local variables to aggregate func- 
tions. 

Definition 1 Let f(S) be an aggregate function. A 
variable, X, is a local variable to f(S) if and only if 
X appears in S and X does not appear in any aggregate 
function that is outer to f{S) or in the NDLP rule that 
contains f{S). 

Definition <jTJ) specifies that every aggregate function 
f(S) has its own set of local variables. For example, for 
the aggregate atom representation of the double sum- 
mation described in Example ([T]), the variable I is local 



variable to the outer aggregate function sum while the 
variable J is local variable to the inner aggregate func- 
tion sum. 

Definition 2 A global variable is a variable that is not 
a local variable. 

Definition 3 An NDLP program, II, is a set of NDLP 
rules of the form 

oi V a 2 V . . . V Ofc f- Ofc+i,. . . ,a m ,not a m+1 , . . . ,not a n 

(4) 
where ax, 02, . • . , Ofc are atoms and a k +x, . . . ,a n are 
atoms or aggregate atoms. 

Let r be an NDLP rule of the form (gj). We 

use head(r) = ai V «2 V . . . V at and body(r) — 
a k+ i, . . . ,a m ,not a m +x, ■ ■ ■ ,not a n . 

2.2 NDLP Programs Semantics 

Definition 4 The local ground instantiation of a sym- 
bolic set S = {F | C} is the set of all local ground pairs 
of the form {9 (F) \ 9 (C)), where 9 is a substitution 
of every local variable appearing in S to a constant from 
U c - 

Definition 5 Let S be a symbolic set. Then, the ground 
instantiation of S is the local ground instantiation of S , 
then followed by the local ground instantiation of every 
symbolic set, S' , appearing in S, then followed by the 
local ground instantiation of every symbolic set, S" , ap- 
pearing in S' , then etc. 

Example 2 Consider the grounding of the aggregate 
function representation of the double summation pre- 
sented in Example fip and described by ^3j), where n = 
m = 2 and with the following facts added. 

d(l,l). d(l,2). d(2,l). d(2,2). e(3,l). 

Since the aggregate function in j3|) contains two levels of 
nesting of the aggregate function sum, then the ground- 
ing is achieved in two steps. The first step is described 
as follows, where A\, A%, and A3 are variables act as 
place holders that are replaceable by constants. 

sum{ {Ax I (1 < 1 < 2), 

sum{ 1 + J I d(l, J), (1 < J < 2)} = Ax), {A 2 \ (1 < 2 < 2), 

sum{ 2 + J I d(2, J), (1 < J < 2)} = A 2 ), {A 3 \ (1 < 3 < 2), 

sum{3 + J I d(3, J),(l < J<2)} = A 3 ) } 

followed by the second step which is described as: 

sum { { Ax I (1 < 1 < 2), 

sum { (1 + l|d(l, 1), (1 < 1 < 2)), (1 + 2|d(l, 2), (1 < 2 < 2)), 

(l + 3|d(l,3),(l<3<2))} = ^ 1 ), 

( A 2 I (1 < 2 < 2), 
sum { (2 + l|d(2, 1), (1 < 1 < 2)), (2 + 2|d(2, 2), (1 < 2 < 2)), 
(2 + 3|d(2,3),(l<3<2)}} = 4 2 ), 

( A 3 I (1 < 3 < 2), 
sum { (3 + l|d(3, 1), (1 < 1 < 2)), (3 + 2|d(3, 2), (1 < 2 < 2)), 
(3 + 3|d(3,3),(l<3<2))} = A 3 ) 

} 



Definition 6 A ground instantiation of an NDLP rule, 
r, is the replacement of each global variable appearing 
in r to a constant from Uc, then followed by the ground 
instantiation of every symbolic set, S , appearing in r. 

The ground instantiation of an NDLP program, LT, is the 
set of all possible ground instantiations of every NDLP 
rule, r, in LT. 

Let X be a set of objects. Then, we use 2 to de- 
note the set of all multisets over elements in X. The 
semantics of the aggregate functions min, max, count, 
sum, and times are defined by the mappings; min, max : 

(2 K - 0) -)• R; count : 2 Uc -> N, sum : 2 R -> E, and 

times : 2 — »■ R, where R is the set of all real num- 
bers, N is the set of all natural numbers, and Uc is the 
Herbrand universe. The application of sum and times 
on the empty multiset returns zero and one respectively. 
The application of count on the empty multiset returns 
zero. However, the application of max and, min on the 
empty multiset is undefined. Let 1 be a symbol that 
does not occur in any NDLP program. 



Definition 7 An interpretation is a subset of the Her- 
brand base Be- 

An atom, a, is true (satisfied) with respect to an inter- 
pretation, /, if a belongs to /; but it is false (unsatisfied) 
otherwise. The negation of an atom, not a, is true (sat- 
isfied) with respect to / if a does not belong to /; but it 
is false (unsatisfied) otherwise. Similarly, the evaluation 
of an aggregate function, and hence the truth valuation 
of an aggregate atom, are established with respect to 
a given interpretation, /, as described by the following 
definitions. 

Definition 8 Let f(S) be a ground aggregate function 
and I be an interpretation. Then, we define Si to be 
the multiset constructed from elements in S , where Si = 
{{F G | (F G | C G ) e SA C G is true in I}. 

Definition 9 Let f(S) be a ground aggregate function 
and I be an interpretation. Then, the evaluation of f(S) 
with respect to I is, f(Si), the result of the application 
of f to Si. f{Si) — J_ if Si is not in the domain of f. 

Example 3 Let L = {d(l, l),d(l,2),d(2, l),d(2,2),e(3, 1)} 

be an interpretation. Thus, the evaluation of the ground 
aggregate function representing the double summation 
in Example f5|) is evaluated w.r.t. L as follows in two 
steps. Considering only the relevant possible values of 
the variables A\, A%, and A3 in Example H|), the first 
step is to evaluate the inner instance of the aggregate 
function sum as follows: 

sum {( 5 I (1 < 1 < 2), sum { 2, 3 } = 5 ), ( 7 | (1 < 2 < 2), 
sum{ 3,4} = 7)} 

followed by evaluating the second instance of the aggre- 
gate function sum as: 

«wi{5,7} = 12 = X£=iE?=i(*+i) 



Definition 10 Let f(S) ~< T be a ground aggregate 
atom and I be an interpretation. Then, f(S) -< T is 
true (satisfied) with respect to I if and only if f(Si) =/= _L 
and f(Sr) -< T. Furthermore, not f(S) -< T is true 
(satisfied) with respect to L if and only if f(Si) = _L or 
f(Si) ± 1 and f(Si) -AT. 

Definition 11 Let II be a ground NDLP program, r be 
a ground NDLP rule of the form O), and L be an inter- 
pretation. Then, 

• L satisfies head(r) iff 3i (1 < i < k) such that L 
satisfies a^. 

• L satisfies body(r) iff\/(k + 1 < i < m) L satisfies 
di and V(m + 1 < j < n) I satisfies not aj . 

• / satisfies r iff I satisfies head(r) whenever L satis- 
fies body(r) or L does not satisfy body(r). 

• L satisfies LT iff L satisfies every NDLP rule, r, in 

n. 

2.3 Answer Sets 

A model for an NDLP program, LT, is an interpretation 
that satisfies LT. A model / of LT is (--minimal if and 
only if there does not exist a model I' of LT such that 

I'd. 

Definition 12 Let LT be a ground NDLP program, r be 
an NDLP rule in LT, and L be an interpretation. Let 
L \= bodyOr) denotes that L satisfies body(r). Then, the 
reduct, II , of II w.r.t. L is the ground NDLP program 
LT 7 where 

n J = {head(r) 4- body{r) \ r £ U A L (= body{r)} 

The reduct H 1 of H w.r.t. / excludes all rules r G LT 
whose body, body(r), is not satisfied by /. The satis- 
faction of body(r) in the definition of the reduct does 
not distinguish between atoms or aggregate atoms or the 
negation of atoms or the negation aggregate atoms. This 
means that if dissatisfaction of body(r) is due to unsatis- 
fied atom or aggregate atom or unsatisfied negated atom 
or negated aggregate atom the consequence is the same, 
which is the exclusion of r from the reduct of LT. 

Definition 13 An interpretation, L , for a ground 
NDLP program, LT, is an answer set for LT if I is C- 
minimal model for LT . 

Observe that the definitions of the reduct and the an- 
swer sets semantics for NDLP programs are general- 
izations of the definitions of the reduct and the an- 
swer sets semantics for a single aggregate over a vari- 
able disjunctive lo gic programs, DLP- 4 , described in 
|Faber et al., 2010] , and hence, generalizations of the 
definitions of the reduct and the answer sets semantics 
for the o riginal disjunctive logic prog rams, DLP, pre- 
sented in [Gelfond and Lifschitz, 1991 . 



2.4 Semantics Properties 

In this section we study the semantics properties of 
NDLP and its relationship to the answer set semantics of 



a single variable a nd a single aggrega te disjunctive logic 
programs, DLP- 4 |Faber et al, 2010] , and its relation- 
ship to the answer set semantics of the original disjunc- 
tive logic programs, DLP Gelfond and Lifschitz, 1991 . 



Theorem 1 Let II be an NDLP program. The answer 
sets of II are Q-minimal models for II. 

The following theorem shows that the answer set se- 
mantics of NDLP subsumes the answer set semantics of 
DLP- 4 



Faber et al., 2010 



and consequently subsumes 
the original answer set semantics of the original disjunc- 
tive logic programs DLP [Gelfond 'a nd Lifschitz, 1991 . 
DLP" 4 programs are NDLP programs with a single ag- 
gregation over a single variable which, unlike NDLP pro- 
grams, do not allow aggregations over function terms. 
DLP programs are NDLP programs without any aggre- 
gate atoms. 

Theorem 2 Let II be a DLP program and L be an in- 
terpretation. Then, I is an answer set for H iff I is an 
answer set for II according to the answer set semantics 
of \Faber et al, 201 Of . 

Proposition 1 Let H be a DLP program and I be an 
interpretation. Then, I is an answer set for H iff I is an 
answer set for II according to the answer set semantics 
of \Gelfond and Lifschitz, 199lf . 

3 Nested Aggregates Answer Set 
Optimization 

In this section, we introduce the syntax and semantics of 
the answer set optimization programs with preferences 
that involve nested multiple aggregates, called nested ag- 
gregates preferences, and denoted by NASO programs, 
that modify and generalize the syntax and semantics of 
the answer set optimization with aggregate preferences 
presented in |Saad and Brewka, 2011 , from a single ag- 
gregate over a single variable preferences to nested mul- 
tiple aggregates over multiple variables and nested mul- 
tiple aggregates over functions that involve multiple vari- 
ables preferences. An NASO program is a logic program 
under the answer set semantics whose answer sets are 
ranked according to preference relations represented in 
the program. 

An NASO program, II, is a union of two sets of 
logic rules II = R gen U R pre f- The first set of logic 
rules, Rgen, is called the generator rules that gener- 
ate the answer sets that satisfy every rule in R gen - 
Rgen is any set of logic rules with well-defined an- 
swer set semantics includ ing normal, extended, and 
disjunctive sets of rules IGelfond and Lifschitz, 1988 



|Gelfond and Lifschitz, 1991] " Faber et al., 2010 , as well 
as disjunctive logic with nested multiple aggregates sets 
of rules presented in the first part of this paper. The sec- 
ond set of logic rules, R pre f, is a set of logic rules that 
represent the user preferences over the answer sets gener- 
ated by Rgen , called the preference rules. The preference 
rules in R pre f are used to rank the generated answer sets 
from the most preferred answer set to the least preferred 



one. An advantage of NASO is that R gen and R pre f 
are independent. This makes preference elicitation eas- 
ier and the whole approach is more intuitive and easy to 
use in practice. We focus on the syntax and semantics of 
the preference rules R pre f of the NASO programs LT = 
Rgen U R P ref , since the syntax and semantics of R gen is 
the same as syntax and semantics of any set of logic rules 
with answer set semantics Gelfond and Lifschitz, 1991[ 



IGelfond and Lifschitz, 1988||Faber et al, 2010| . 



3.1 NASO Programs Syntax 

The language of NASO programs is the same as the lan- 
guage NDLP programs, presented in the first part of 
this paper, except that in the language of NASO pro- 
grams classical negation is allowed. Let £ be a first-order 
language with finitely many predicate symbols, function 
symbols, constants, and infinitely many variables. A 
term is a constant, a variable or a function. A literal 
is either an atom a in Be or the negation of a (~<a), 
where Be is the Herbrand base of C and -> is the classi- 
cal negation. The Herbrand universe of C is denoted by 
V~c- Non-monotonic negation or the negation as failure 
is denoted by not. Let Lit be the set of all literals in C, 
where Lit = {a\a £ Be} U {^a|a £ Be}- 

An expression of the form {F | C} is called a sym- 
bolic set, where F is a single variable or a function 
F = g(Xi, . . . , X n ), where X\, . . . , X n are variables or 
functions terms and C is a conjunction of literals, non- 
monotonic negation of literals, aggregate atoms, and 
non-monotonic negation of aggregate atoms (defined be- 
low). However, a ground set is a set of pairs of the 
form (Fq I Cg) such that Fq is a constant and Cq is a 
conjunction of ground literals, non-monotonic negation 
of ground literals, ground aggregate atoms, and non- 
monotonic negotiation of ground aggregate atoms (de- 
fined below). A ground set or a symbolic set is called a 
set term. We say f(S) is an aggregate function if / is 
an aggregate function symbol and S is a set term, where 
/ £ {min, max, count, sum, times}. If f(S) is an aggre- 
gate function and T is a constant, variable or function 
term called guard, then we say f(S) ~< T is an aggre- 
gate atom, where ^<G {=, A <,>,<,>}• An optimiza- 
tion aggregate is an expression of the form max(f(S)) 
or min(f(S)), where S is a set term and / is an ag- 
gregate function symbol. Let A be a set of literals, ag- 
gregate atoms, and optimization aggregates. A boolean 
combination over A is a boolean formula over literals, 
aggregates atoms, and optimization aggregates in A con- 
structed by conjunction, disjunction, and non-monotonic 
negation (not), where non-monotonic negation is com- 
bined only with literals and aggregate atoms. 

Similar to NDLP, an aggregate function f(S) in an 
aggregate atom f(S) -< T can be defined in terms of 
other aggregate atoms f'(S') ~< T" and hence in terms 
of other aggregate functions f'(S'). This is because the 
conjunction in the definition of the symbolic set S of the 
aggregate function f(S) in the aggregate atom f(S) -< T 
can contain other aggregate atoms as constituents. This 
implies that multiple level of nested aggregate functions 



calculation is achieved by allowing aggregate functions 
of the aggregate atoms to be defined in terms of other 
aggregate atoms. 

Let f(S) be an aggregate function. A variable, X, is 
a local variable to f(S) if and only if X appears in S 
and X does not appear in any aggregate function that 
is outer to f(S) or in the preference rule the contains 
f(S). A global variable is a variable that is not a local 
variable. 

Definition 14 A preference rule, r, over a set of liter- 
als, aggregate atoms, and optimization aggregates, A, is 
an expression of the form 

C\ >- C 2 y . . . y Ck 4- h+i, • • • , l m , not l m +i, ■ ■ ■ , not l n 

(5) 
where lk+i, ■ • ■ ,l n are literals or aggregate atoms and 
C\, C 2 , . . • , Ck are boolean combinations over A. 

Let r be a preference rule, body(r) = 
lk+i,...,l m ,not l m+ i,...,nat l n , and head(r) = 
C\ >- C 2 >~ ••• >~ Ck- Intuitively, a preference rule, r, 
says that any answer set that satisfies body(r) and C\ 
is preferred over answer sets that satisfy body(r), some 
Ci (2 < i < k), but not C\, and any answer set that 
satisfies body(r) and C'2 is preferred over answer sets 
that satisfy body, some Ci (3 < i < k), but neither C\ 
nor C 2 , etc. 

Recalling, an NASO program is a union of two sets 
of logic rules II = R gen U R pre f, where R gen is a set of 
logic rules with answer set semantics, called the genera- 
tor rules, and R pre f is a set of preference rules. 

3.2 NASO Programs Semantics 

In defining the semantics of aggregate functions and ag- 
gregate atoms for NASO programs as they syntactically 
defined in the previous section, we follow the same se- 
mantics of aggregate functions and aggregate atoms as 
they defined for NDLP programs. We use for NASO 
programs the same notions of local ground instantiation 
and ground instantiation of symbolic sets as they defined 
for NDLP programs. Similarly, the semantics of the ag- 
gregate functions min, max, count, sum, and times are 
defined by the same mappings as defined for NDLP pro- 
grams. Let i be a symbol that does not occur in any 
NASO program. 

A ground instantiation of a preference rule, r, is the 
replacement of each global variable appearing in r to a 
constant from Uc, then followed by the ground instantia- 
tion of every symbolic set, S, appearing in r. The ground 
instantiation of an NASO program, II = R gen L)R pre f, is 
the set of all possible ground instantiations of every rule, 
r, in II. Let I be an answer set for R gen in an NASO 
program, II = R gen U R pre f, f(S) -<Tbca ground ag- 
gregate atom, and Si = fF G \ {F G \ Cq) € S and Cq 
is true in /]}• be the multiset constructed from elements 
in 5. Then, the evaluation of f(S) with respect to the 
answer set 7 is f(Si), where f(Si) = _L if Si is not in 
the domain of /. 



Definition 15 Let H = R gen UR pre f be a ground NASO 
program, I be an answer set for R gen , and r be a pref- 
erence rule in R pre f ■ Then the satisfaction of a boolean 
combination, C, appearing in head(r), by I, denoted by 
I \= C , is defined inductively as follows: 

• I^=liffl€l. 

• I \= not I iff I <£ I. 

. I h f(S) -< T iff f (Si) ± J_ and f(Si) -< T. 

• 7 |= not f(S) <T iff f(Si) = _L or f(Si) ^ J_ and 
f(Si) it T. 

• 7 |= max(f(S)) iff f (Si) ^ _L, and for any answer 
set I', f(S r ) ± J_ andf(S r ) < f(Sj) or f(Si) ± JL 
and f(Si') = _L. 

• 7 |= min(f(S)) iff f (Si) 7^ _L ; and for any answer 
set I', f(S r ) £ J_ and f (Si) < f(S r ) or f(S T ) ± JL 
and f(Si') = _L. 

• I [= d A C 2 iff I \= Ci and I f= C 2 . 

• /HCiVC 2 iffl^Ci orI^C 2 . 

The application of any aggregate function, /, except 
count, on a singleton {a;}, returns x, i.e., f({x}) = x. 
Therefore, we use max(S) and min(S) as abbrevia- 
tions for the optimization aggregates max(f(S)) and 
min(f(S)) respectively, where S is a singleton and / 
is any arbitrary aggregate function except count. 

Definition 16 Let H = R gen LlR pre f be a ground NASO 
program, 7 be an answer set for R gen , and r be a prefer- 
ence rule in R pre f ■ Then the satisfaction of the body of 
r by I, denoted by I \= body(r), is defined inductively as 
follows: 

• I\=liffl€l. 

• I \= not I iff I <£ I. 

. 7 h f(S) -< T iff f (Si) ^ J_ and f(S T ) ■< T. 

• I\=not f(S) ~< T iff f(Si) = ±or f(Si) ± J_ and 
f(Si) ^ T. 

• 7 |= body(r) iff V(fc + 1 < i < m), 7 |= li, and 
V(m + 1 < j < n), L \= not lj . 

The following definition specifics the satisfaction of the 
preference rules. 

Definition 17 LetH = R gen L)R pre f be a ground NASO 
program, I be an answer set for R gen , r be a preference 
rule in R pre f , and Ci be in head(r). Then, we define the 
following notions of satisfaction of r by I: 

• I |=j r iff 7 |= body(r) and I \= d. 

• 7 \=i rr r iff L \= body(r) and I does not satisfy any 
of Ci in head(r) . 

• 7 \=i rr r iff I does not satisfy body(r). 

Definition 18 LetH = R gen L)R pre f be a ground NASO 
program, 7i,7 2 be two answer sets of R gen , r be a pref- 
erence rule in R pre f , and Ci be boolean combination ap- 
pearing in head(r). Then, L\ is strictly preferred over 
7 2 w.r.t. r, denoted by I\ y r 7 2; iff one of the following 
holds: 



• I\ |=j r and Ii \=j r and i < j , 

where i — min-{7 | I\ |=/ r} and j ~ Tam{l | I2 \=i 
r}. 

• I 1 \= i r and I 2 \=i rr r. 

We say, I\ and 1 2 are equally preferred w.r.t. r, denoted 
by I\ = r I2, iff one of the following holds: 

• I\ \=i r and I2 |=j r, where i = mhx{l \ I\ \=i r} = 
min{l I I 2 \=i r }- 

• h \=irr r and I 2 \=i rr r. 

We say, I\ is at least as preferred as I2 w.r.t. r, denoted 
by h hr h, iff h >r h or h =,. I 2 - 

Definition (fT5|) specifies the ranking of the answer sets 
according to a preference rule. The following definitions 
characterize the ranking of the answer sets with respect 
to a set of preference rules. 

Definition 19 (Pareto Preference) Let II = R gen U 
Rpref be an NASO program and I\,l2 be answer sets 
of Rgen- Then, I\ is (Pareto) preferred over I2 w.r.t. 
Rpref , denoted by I\ ^n . I2, iff there exists at least 
one preference rule r 6 R pre f such that I\ >- r I2 and for 
every other rule r' € R pre f, I\ ^>' ^2- We say, I\ and 
I2 are equally (Pareto) preferred w.r.t. R pre f, denoted 
by h =R pre/ h, iff for all r € Rpref, h =r h- 

Definition 20 (Maximal Preference) Let 

II = Rgen U Rpref be an NASO program and I\,I% 
be answer sets of R gen - Then, L\ is (Maximal) preferred 
over L2 w.r.t. R pre f , denoted by L\ >-_r p ,. c/ I2, iff 

\{r € Rpref\h h r h}\ > \{r € Rpref\h hr h}\- 

We say, I\ and 1% are equally (Maximal) preferred w.r.t. 
R pre f, denoted by h = flpre/ I 2 , iff 

\{r € Rpref \h hrh}\ = \{re Rpref \h tr h}\- 

Observe that the Maximal preference relation is more 
general than the Pareto preference relation, since the 
Maximal preference definition subsumes the Pareto pref- 
erence relation. Under the Pareto preference rela- 
tion, the following result shows that the syntax and 
semantics of NASO programs subsume the syntax 
and semantics of the answer set optimization pro- 

since there is no no- 



grams of Brewka et al., 2003 



tion of Maximal preference relation was introduced in 



Brewka et al., 2003 . This is assuming that the answer 
set optimization programs of Brewka et al., 2003 assign 



the lowest rank to the answer sets that do not satisfy nei- 
ther the body nor the head of preference rules. 

Theorem 3 Let II = R gen U Rpref be an NASO pro- 
gram without either aggregate atoms or optimization 
aggregates and I\ , I2 be answer sets of R 



gen- 



The 



L\ is Pareto preferred over L2 w.r.t. R pre f iff Ii is 
Pareto preferred over L2 w.r.t. Rpref according to 
\Brewka et al., 2003J . 

However, the following result shows that the syn- 
tax and semantics of NASO programs subsume the 



syntax and semantics of a single aggregate over 
a single variable answer set optimization programs 
of |Saad and Brewka, 2011] , under both the Pareto 
and the Maximal preference relations, since the no- 
tion of Maximal preference relation introduced in 
|Saad and Brewka, 201 1| is a special case of the Maxi- 
mal preference relation presented in this paper. 

Theorem 4 Let H — R gen U Rpref be a single aggre- 
gate over a single variable NASO program and ii,/2 
be answer sets of R gen - Then, I\ is Pareto (Maxi- 
mal) preferred over I2 w.r.t. R pre f iff I\ is Pareto 
(Maximal) preferred over L2 w.r.t. R pre f according to 
\Saad and Brewka, 20 ll 1 . 



4 Probabilistic Traveling Salesman 
Problem 

In this section we show that any instance of a priori 
optimization probabilistic traveling salesman problem 
(PTSP) |Jaillet, 1988| can be intuitively and easily rep- 
resented and solved by the framework of nested multiple 
aggregates answer set optimization programs. In addi- 
tion, we show that with a minor modification to our rep- 
resentation of any instance of PTSP problem in NASO, 
we can intuitively solve a corresponding classical travel- 
ing salesman problem (TSP) in our framework in par- 
ticular, and in answer set programming in general, since 
finding the optimal tour for TSPs has not been consid- 
ered before in answer set programming literature. This 
is because the emphasis in solving TSPs in answer set 
programming literature was on generating the possible 
tours of a given TSP rather than finding the optimal 
tour for that TSP, which is the tour with the minimum 
length. The reason for that is tours of a TSP are rep- 
resented in answer set programming as answer sets and 
answer set programming is incapable of reasoning across 
answer sets to find the answer set that represents the 
tour with the minimum length. Therefore, a different 
framework that is capable of reasoning across answer 
sets is required, this framework is NASO programs. 

Observe that our representation of PTSPs and conse- 
quently of TSPs is built on top of the existing answer 
set programming representation of TSPs, whose aim is 
to find all possible tours. This makes our NASO pro- 
gram representation further intuitive and keeps inline 
with the existing body of work in answer set program- 
ming. In addition, it shows and gives insight that many 
optimization problems can be solved by NASO frame- 
work in the same way by similar intuitive modifications 
to the existing answer set programming representation 
of corresponding similar problems. 

4.1 Probabilistic Traveling Salesman 
Problem in NASO 

Recalling, a Probabilistic Traveling Salesman Problem 
(PTSP), G, is a tuple of the form G = (V,E,D,P), 
where V is a set of n vertices, E is a set of edges, D 
is an n x n distance matrix, where dy is distance from 



a vertex j to a vertex j, and P = {pi,ps, ■ ■ ■ ,Pn} is a 
probability distribution over vertices in V, where pi is 
the probability a vertex i must be visited. The aim is to 
find a priori a tour, t — (0, 1, 2, 3, 4, ... , n, 0), of G, with 
the minimum expected length, where is the starting 
point from which the tour always starts and ends and 
whose probability po = 1. The expected length, E[L t ], 
of a tour t is given by formula (JXJ) . 

Any Probabilistic Traveling Salesman Problem 
(PTSP), G = (V,E,D,P), is represented as an NASO 
program II = R gen U R pre f, where the generator rules 
in R gen generate the answer sets the represents all the 
possible tours in G and the preference rules in R pre f 
rank the answer sets that correspond to all the possible 
tours in G from the tour with the minimum expected 
length to the tour with the maximum expected length. 
The representation proceeds as follows. 

• Every vertex x in V whose probability is p in P is 



represented in R gen as a fact of the form 
vertex{x,p) <— 



(G) 



In addition to the starting vertex y € V is repre- 
sented in Rg en as a fact of the form 

start(y) <— 

Every edge (x, y) in E is represented in R gen 
fact of the form 



edge(x,y) 



(7) 



(8) 



Every element d xy in the distance matrix D is rep- 



resented in Rgen as a fact of the form 



distance(x,y,d xy ) <— 



(9) 



The logic rules that generate all the possible tours 
from, G, are represented in R gen by the disjunctive 
logic rules: 

rinTour(X, Y) V outTour(X, Y) <- start(X), 

edge(X,Y). (10) 

inTour{X, Y) V outT our(X , Y) ^— reached(X) , 

edge(X,Y). (11) 

reached(Y) <- inTour(X,Y). (12) 

«- inTour(X,Y),inTour(X,Y 1 ),Y ^ Y x . (13) 

<- inTour(X, Y),inTour(X 1 ,Y), X ^ X x . (14) 

<— vertex(X,P),not reached(X). (15) 

(16) 

The indexing of vertices in any given tour, t, by their 
order of appearance in the tour t is represented in 
Rgen by the logic rules: 

index(0,X) <- start(X). (17) 
index(I + 1, Y) «- inTour(X, Y), edge(X, Y), 

index(I,X). (18) 



• The computation of the expected length of each tour 
and the ranking of the answer sets corresponding to 
the tours in G from the minimum expected length 
tour to the maximum expected length tour is repre- 
sented in R pr ef by the preference rule (denoted by 
the preference rule opt): 

min( sum { A 
index(I,X), vertex{X,Px), < I < n, 
sum { Px * Py * D * M 
index(J,Y), vertex(Y, Py), distance{X,Y,D), 
7 + 1 < J<n + 1, 
times { (1 — Pz) 
index(K,Z), vertex(Z,P z ), I + 1 < K < J - 1} = M } 

)<- 
Observe that the logic rules (fTU]) - ([T2l are exactly a typ- 
ical set of logic rule s that are used in answer set pro- 
gramming literature |Eiter et at, 2000] to find Hamilto- 
nian cycles of a given graph adapted to deal with graphs 
with probabilistic vertices. It can be easily seen that the 
nested aggregate function within the optimization ag- 
gregate min in the above preference rule is exactly the 
representation of the expected length of a tour described 
by formula (J]). 

Theorem 5 Let G = (V, E, D, P) be a Probabilistic 
Traveling Salesman Problem and II = R gen U R pre f be 
the NASO program representation of G. Then, a tour, 
t, is the optimal tour for G with the minimum expected 
length iff an answer set, I, of R gen is the top preferred 
answer set w.r.t. R pre f. 

4.2 Traveling Salesman Problem in NASO 

In this section, we show that with a minor modification 
to the NASO program representation of a probabilis- 
tic traveling salesman problem, we can intuitively solve 
a corresponding classical traveling salesman problem 
(TSP). A classical Traveling Salesman Problem (TSP), 
G, is a tuple of the form G = (V,E,D), where V is a 
set of n vertices, E is a set of edges, and Disannxn 
distance matrix, where d(i,j) is an element in D that 
represents the distance from a vertex i to a vertex j. 
The aim is to find a tour, t — (0, 1,2,3,4,..., n, 0), of 
G, with the minimum length, where is the starting 
point from which the tour always starts and ends and 
the length, L t , of a tour t is calculated by 



L t =^2d(i,i + l) 



(19) 



where n + 1 = 0. A Traveling Salesman Problem (TSP), 
G = (V,E,D), is represented as an NASO program, 
II = R gen U Rpref , where the generator rules in R gen gen- 
erate the answer sets the represents all the possible tours 
in G and the preference rules in R pre f rank the answer 
sets that correspond to all the possible tours in G from 
the tour with the minimum length to the tour with the 
maximum length. The generator rules in R gen consists 
of the logic rules from the logic rule ([6]) through the logic 
rule (|18p after replacing vertex(x,p) in the logic rule (J5]) 



= ^4} 



by vertex(x) and replacing v ertex(X 1 P) in the logic rule 
(fT5| by vertex(X). The computation of the length of 
each tour and the ranking of the answer sets correspond- 
ing to the tours in G from the minimum length tour to 
the maximum length tour is represented in R pre f by the 
preference rule: 

min( sum { D \ distance(X,Y,D), index(I,X), 
index(I + 1,Y), < I < n } ) <r- 

Theorem 6 Let G = {V, E, D) be a classical Traveling 
Salesman Problem and IT = R gen U R pre f be the NASO 
program representation of G. Then, a tour, t, is the 
optimal tour for G with the minimum length iff an an- 
swer set, I , of R gen is the top preferred answer set w.r.t. 

ftpref • 

4.3 Probabilistic Traveling Salesman 
Problem in NDLP 

We also show that a Probabilistic Traveling Salesman 
Problem, G = (V,E,D,P), can be represented as an 
NDLP program, II, whose answer sets corresponds to 
tours in G. Although, NDLP programs framework is still 
capable of finding all possible tours in G, represented as 
answer sets, and computing the expected length of each 
tour, the semantics of NDLP programs is not able to lo- 
cate the answer set that corresponds to the optimal tour 
of G, since it does not have the capability of reasoning 
across the answer sets. However, PTSPs can be repre- 
sented and solved in NDLP programs framework in two 
steps. The first step is to represent a probabilistic trav- 
eling salesman problem, G, by an NDLP program whose 
answer sets correspond to all possible tours in G along 
with their expected length. The second step is to deter- 
mine the minimum expected length optimal tour, t, in G 
represented by an answer set of the NDLP program rep- 
resentation of G by means of any appropriate procedure 
internal or external to NDLP programs framework. 

The NDLP program, II, representation of a Proba- 
bilistic Travelling Salesman Problem, G, consists of the 
NDLP rules from the logic rule ^ through the logic rule 
(JT5J), in addition to the following NDLP rule, where X in 
length(X) in the NDLP rule below represents the value 
of the expected length of a tour. 

length(X) <— sum { A \ 

index(I,X), vertex{X,Px), < I < n, 

sum { Px * Py * D * M \ 

index(J,Y), vertex(Y, Py), distance(X,Y, D), 

I+l<J<n+l, times { (1 - P z ) \ 

index(K, Z), vertex(Z, P z ), I + 1 < K < J - 1} = M } 

= A} = X. 

Theorem 7 Let G = (V, E, D, P) be a Probabilistic 
Traveling Salesman Problem and II be the NDLP pro- 
gram representation of G. Then, t is a tour for G iff L 
is an answer set for II corresponds to t, where the ex- 
pected length, E[L t ], oft is equal to the value of X in 
length(X) that is satisfied by I . 



One possible way to find the optimal tour with the 
minimum expected length is to add the following two 
NDLP rules to the NDLP program representation, II, of 
a probabilistic traveling salesman problem problem, G, 
where X in min(X) in the NDLP rules below represents 
the minimum expected length value of a tour 

min(X) <— length(X) , not maximal(X). (20) 
maximal(X) 4— length(X),length(Xi), X\ < X. (21) 

Theorem 8 Let G = {V, E, D, P) be a Probabilistic 
Traveling Salesman Problem and II be the NDLP pro- 
gram representation of G. Then, a tour, t, is the optimal 
tour for G with the minimum expected length iff an an- 
swer set L of II is the only answer set of II that satisfies 
min(X). 

However, the addition of the NDLP rules (HOJ) and (12TJ) 
to the NDLP program representation of a probabilistic 
traveling salesman problem leads when grounded to large 
number of ground rules, especially with the graphs that 
involve sufficiently large number of vertices and large 
number of edges, which is always the case with PTSPs. 
Therefore, a simple linear search for the minimum value 
of X in length(X) over all the answer sets of the NDLP 
program representation of a probabilistic traveling sales- 
man problem is likely to be more efficient than perform- 
ing the same task using logic rules. 

4.4 Working Example 

This section shows that a probabilistic traveling sales- 
man problem instance can be intuitively and easily rep- 
resented and solved by the NASO programs framework. 
Consider this instance of the probabilistic traveling sales- 
man problem, G = (V,E,D,P), where (V,E) is com- 
plete undirected graph where V — {a, b, c, d} and E = 
{{a, b), (a, c), (a, d), (6, c), (b, d), (c, d), (b, a), (c, a), (d, a), 
(c,b),(d,b),(d,c)}, the probability distribution, P, is 
given as p a — 1, pb — 0.3, p c = 0.7, and pd — 0.4, and 
the distance matrix, D, is given by d a b = 40, d ac = 40, 
d ad = 22, d bc = 40, d bd = 25, d cd = 22, d ba = 4o| 
d ca = 40, d da = 22, d cb = 40, d db = 25, and d dc = 22. 

Let n = R gen \J Rp re f be the NASO program represen- 
tation of this probabilistic traveling salesman problem 
instance, G, where R pre f consists of the preference rule 
opt, with n = 3, and R gen consists of the logic rules (I10[) 
through HH), in addition to the fats: 

vertex(a,l). vertex(b,0.3). vertex(c, 0.7). vertex(d, 0.4). 

edge(a,b). edge(a,c). edge(a,d). edge(b,c). edge(b,d). 

edge(c,d). edge(b,a). edge(c,a). edge(d,a). edge(c,b). 

edge(d,b). edge(d,c). distance(a,c,40). distance(a,b,40). 

distance(a,d, 22). distance(b,c,AQ). distance(b,d,25). 

distance(c,d, 22). distance(c, a, 40). distance(b,a,A0). 

distance(c,b,40). distance(d,b,25). distance(d, a, 22). 

distance(d,c, 22). start(a). 

By considering only the relevant atoms, R gen , has six 
answer sets which correspond to the six available tours 
in G These answer sets are: 
I\ = {start(a), inTour(a,d), inTour(d,b), inTour(b,c), 
inTour(c,a), index(0,a), index{l,d), index(2,b), 
index(3,c), index(4,a)} 



Ii = {start(a), inTour(a,d), inTour(d,c), inTour{c,b), 
inTour{b,a), index(0,a), index{\,d), index(2,c), 
index{3,b), index(A,a)} 

13 = {start(a), inTour(a,b), inTour(b,d), inTour(d,c), 
inTour(c,a), index(0,a), index(l,b), index(2,d), 
index(3,c), index(A,a)} 

14 = {start(a), inTour{a,c), inTour(c,b), inTour(b,d), 
inTour(d,a), index(0,a), index(l,c), index(2,b), 
index(3,d), index(4,a)} 

15 = {start(a), inTour{a,c), inTour{c,d), inTour{d,b), 
inTour{b,a), index(0,a), index(l,c), index(2,d), 
index(3,b), index(A,a)} 

Iq = {start(a), inTour(a,b), inTour(b,c), inTour{c,d), 
inTour{d,a), index(0,a), index(l,b), index(2,c), 
index(3,d), index(4,a)} 

There are six tours for, G, which are ti, t%, £3, £4, £5, 
and te that correspond to the six answer sets I\, I2, I3, 
I4, I5, and Iq respectively. The expected length, E[L ti ], 
of each tour, t%, in G, calculated by formula ([1]), is given 
as: 

E[L tl ] = 75.16, E[L t2 ] = 76.67, E[L t3 ] = 76.92, 
E[L U ] = 76.92, E[L t5 ] = 76.92, E[L te } = 76.67 

that exactly corresponds to the evaluation of the nested 
aggregate function within the optimization aggregate, 
min, of the preference rule opt contained in the set of 
the preference rules R pre f of the NASO program repre- 
sentation, II = R gen U Rp re f , of the probabilistic travel- 
ing salesman problem instance G. It is clear that t\ is 
the tour with minimum expected length, E[L tl ] — 75.16, 
and hence the optimal tour for G, as well as, I\ is the 
top answer set of R gen with respect to R pre f. This is 
because 

h Hi °V^i h \=irr Opt, I 3 ^ irr Opt, I4 ^ irr Opt, 
h \=irr Opt, I 6 ^ irr Opt. 

Similarly, this probabilistic traveling salesman problem 
instance can be represented by an NDLP program, II', 
that consists of the NDLP rules (JTOJl through ([2"T]). as 
well as the NDLP rule representation of the expected 
length (the definition of the length(X) predicate), and 
by considering only the relevant atoms, it can be easily 
seen that 7i is the only answer set of LT' that satisfies 
min(X) and coincides with the optimal minimum ex- 
pected length tour of G. 

5 Conclusions and Related Work 

We developed a generalization for the syntax and seman- 
tics of both answer set optimization programming and 
answer set programming from a single aggregate over a 
single variable to nested multiple aggregates over mul- 
tiple variables and to nested multiple aggregates over 
functions involving multiple variables. In addition, we 
showed the applicability of both the nested aggregates 
answer set optimization programming and the nested 
aggregates answer set programming to the a priori op- 
timization Probabilistic Traveling Salesman Problem, a 



fundamental stochastic optimization problem in Oper- 
ation Research. Moreover, we showed that with mi- 
nor modification to the probabilistic traveling salesman 
problems representation in nested aggregate answer set 
optimization programming and answer set optimization 
programming, a corresponding classical traveling sales- 
man problems can be intuitively represented and solved 
in both frameworks. 

We showed that the nested aggregates answer set op- 
timization programming framework presented in this 
paper modifies and generalizes a single aggregate over 
a single variable answer set optimization programming 
framework of [Saad and Brewka, 2011 as well as a gen- 
eralization of the answer set optimization programming 
framework of [Brewka et al., 2003] . We showed that the 
nested aggregates answer set programming framework 
presented in this paper subsumes a single aggregate 
over a single variable answer set programming frame- 
work of Faber et al., 2010 as well as a generalization 
of the original answer set programming framework of 
|Gelfond and Lifschitz, 1991] . 

The major difference in the development presented in 
this paper is that we allow the ability to reasoning in the 
presence of nested multiple aggregates in answer sets in 
general and in both answer set optimization program- 
ming and in answer set programming in particular. The 
only answer set optimization framework that allows ag- 
gregates is Saad and Brewka, 2011] , which is answer set 
optimization programs with a single aggregate over a sin- 
gle variable. The existing answer set programming liter- 
ature considered a single aggregate over a single variable 
answer set programming approaches. A comprehensive 
comparisons of the existing a single aggregate over a sin- 
gle variable answer set programming approaches is found 
in |Faber et al, 2010] . 
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